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Abstract 

Models of a rotating universe have been studied widely since Godel , who showed an example 
that is consistent with General Relativity (GR). By now, the possibility of a rotating universe has 
been discussed comprehensively in the framework of some types of Bianchi's models, such as Type 
V, VII and IX P,!^, and different approaches have been proposed to constrain the rotation. Recent 
discoveries of some non-Gaussian properties of the Cosmic Microwave Background Anisotropies 

(CMBA) |],y,y,Q,y,y,Q, such as the suppression of the quadrupole and the alignment of some 

multipoles draw attention to some Bianchi models with rotation [ill. However, cosmological 
data, such as those of the CMBA, strongly prefer a homogeneous and isotropic model. Therefore, 
it is of interest to discuss the rotation of the universe as a perturbation of the Robertson- Walker 
metric, to constrain the rotating speed by cosmological data and to discuss whether it could 
be the origin of the non-Gaussian properties of the CMBA mentioned above. Here, we derive 
the general form of the metric (up to 2nd-order perturbations) which is compatible with the 
rotation perturbation in a flat A-CDM universe. By comparing the 2nd-order Sachs- Wolfe effect 



m m 



17l | due to rotation with the CMBA data, we constrain the angular speed of 
the rotation to be less than 10^^ rad yr~^ at the last scattering surface. This provides the first 
constraint on the shear-free rotation of a ACDM universe. 
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I. INTRODUCTION 



In contrast to most astronomical objects (such as planets, stars and galaxies), the universe 
as a whole is usually considered to be non-rotational. However, the possibility that the 
universe rotates should not be ignored, since solutions of GR corresponding to a rotating 



universe have been found 



18| indicating that a global rotation is physically allowed. 



Although it is widely assumed that galaxies align randomly and have zero net angular 
momentum, there have been many investigations on the general alignment of galaxies. Some 
even claim that a weak alignment of galaxies does exist 



IS 



used to explain [2 1[| the recently discovered non-Gaussian propert 



Such an alignment n aay be 



the CMBA. Furthermore, Jaffe et al. 
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iesQ, y, fffl y, y, [ 



llQ( of 



12] suggest that the Bianchi Type Vllh model 



with a global rotation may be used to explain some anomalies of the CMBA. In addition, 
;he existence of a global rotation may contradict the inflationary model of the early universe 
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2J] and Mach's Principle 
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Rotational perturbations may also be used to 

Therefore, the study of 



determine whether or not the universe is open or closed [sl, I27I 1 
global rotation is of interest in many different aspects of cosmology, and constraint of the 
rotation speed of the universe is important. 

The most popular approach to constrain the magnitude of the global rotation speed is to 
make use of data from the CMBA because of their precision. Most discussions y,]^, Q, 29] 
have focused on homogeneous cosmological models, i.e. Bianchi models [s^. To be consistent 
with obeservations, only Bianchi Type V, VIIq, Vllh and IX models, which include the 
Robertson- Walker model as a special case, are considered. The constraints of the global 
rotation speed obtained depend on the parameters of the models. Besides, shear and vorticity 
are inseparable in these works [si], i.e. zero shear automatically implies zero vorticity. 

There are many other approaches to constrain the global rotation. Based on the idea 
;hat a global rotation induces a total net spin of galaxies, the global rotation can be limited 
32l |. Moreover, empirical relations between angular momenta and mass of galaxies/clusters. 



such as J ~ 
global rotation 32 



'or spiral galaxies and J ~ for clusters can be explained by the 



33|. The acceleration caused by the global rotation may be used to 



explain parts of the accelerating expansion of our universe, and thus the global rotation 
can be constrained by Supernova Type la data |^3^. Recently, some studies of the CMB 
polarization induced by the global rotation are published 35| providing potential constraints 
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in the future. 

To develop a model that preserves the homogeneity and isotropy of the mean CMB, we 
study the rotation of the universe as a perturbation in the Robert son- Walker framework 
with a cosmological constant in this paper. Unlike the Bianchi models, such an approach 
allows to have non-zero rotation but trivial shear. Since the global rotation does not have any 
influences on the Ist-order Sachs- Wolfe effect (SW effect), we need to calculate the metric up 
to 2nd-order perturbations and the 2nd-order SW effect. Then, we will constrain the angular 
speed of the rotation using recent data on CMBA 36|]. Our model is inhomogeneous with 
an axial symmetry in general. The global rotation in our model is not only time-dependent 
but also radial-dependent. 

II. SOLUTIONS OF THE EINSTEIN FIELD EQUATIONS WITH ROTATIONAL 
PERTURBATIONS 

The line element of a flat rotational universe possesses an axial symmetry and can be 



written in the form of 37 1 



ds^ = g^^dx'^dx" (1) 
= a\i]){[l - f\r,r])]drf - [1 - h{r,r])]dr'^ - [1 - h{r,r])]r^de^ - [1 - k{r,r])]dz^ 
+2r^a{ri)n{r,r])dedr]}, (2) 

where /i and u = 0,1, 2, 3, x'^ = {rj, r, 9, z}, rj is the conformal time defined by dt = a{ri)dr) 
with t the cosmological time, r, 6 and z are the cylindrical coordinates in the comoving 
frame of the universe, z is the axis of rotation, a{ri) is the scale factor of the universe with 
a(?7 = 0) = 1 at the present time, Q{r,ri) is the angular velocity of the metric observed 
from an inertial frame whose origin is on the rotational axis, and /(r, rj), h{r, if) and k{r, t]) 
are the perturbations on the (0, 0), (1, 1), (2, 2), (3, 3)-components of the metric due to 
the rotation. Because of the cylindrical symmetry, the perturbation functions due to the 
rotation are also independent of 9 and z. 

Here, we assume that the norm of g^u—odfiu, where og^^i, is the unperturbed metric, is 
much smaller than that of (7^,^. Explicitly, we assume that the rotation is slow, so that 
ra{'r])Q{r,T]) <^ 1, and we can think of Qrot = maxi7(r, 77) for r, rj within the last scattering 
surface as the perturbation parameter. By parity consideration, we can see that Q{r, rj) is 



composed of only odd powers of Qrot, whereas f{r,ri), h{r,r]) and k{r,r]), being density and 
pressure perturbations, only even powers. Since we are interested only up to second-order 
perturbations, we will consider Q{r,ri) to be first-order and f{r,T]), h{r,T]) and k{r,ri) to be 



second-order. The metric Eq. (1.1) in Ref. 



will be recovered if we truncate ours up to 



the first-order. Since the effect of the rotation on the CMBA is independent of the parity, 
we expect that the SW effect due to rotation occurs in even orders of Qrot only. 

The Einstein Field equations (EFEs) for a universe with cosmological constant A are 

Rf,u - ^QfiuR + ^Qf^u = + ^Qfjiu = SuTf,^, (3) 

where T^,^ = {p+P)u^Uy — Pg^y is the stress-energy tensor for a perfect fluid, R^y is the Ricci 
curvature tensor, R is the scalar curvature, p is the mass-energy density, P is the pressure 
and = a~^{ri)dx^ / drj is the four- velocity of the fluid in the comoving frame. Here, we set 
G = c = l. 

If fi(r, = ^(77), the universe is homogeneous and ^(r/) expresses the angular velocity 
of the universe observed anywhere in the comoving frame. Otherwise, the universe is inho- 
mogeneous and the observer at the rotating axis passing through the origin is distinct. To 
solve Eq. (3) up to 2nd-order in VLrot^ we expand all quantities: 

9^lv = Qg^lu + lg^lu + 2g^ly7 (4) 
P = oP + iP + 2P, (5) 
P = oP + iP + 2P, (6) 

= + lU^ + (7) 

where the subscripts indicate the corresponding orders of perturbations. 
The zeroth-order EFEs give rise to the standard Friedmann equations: 

3d2(r7) 



a?{vi) 



+ Ka\r]) = 87Ta\r])op{v), (8) 



-^a\v) + ^ - ^ = 8na\v)oPiv). (9) 
a^['ri) a[ri) 

Once the equation of state (EOS) of the fluid is given, we can determine the scale factor 
a(?7), the density op('7) and the pressure oPij]) with the equations above. In this paper, we 
consider a universe with oP{v) = 0- However, the following formalism can be applied to any 
fluid with a specified EOS. 



From the temporal-spatial EFEs and the condition Ufj_u^ = 1, we have 0^^(77) = 
{a{ri), 0, 0, 0) and iUo{r, rj) = 0. The first-order EFEs then give: 



ip{r,r]) = 0, iP(r,?7) = 0, iMi(r,?7) = 0, iU3{r,r]) = 0, (10) 
iM2(r,r/) = r^a^{ri)[n{r,ri) - iu^{r,ri)] = ^^^^^^ - 4a{r]yd{r]) l^^^'i^^v) + r^^"ir,ri)]lll) 
= 3d{ri)n'{r,ri) + a{ri)n'{r,r]), (12) 

where the dots refer to derivatives with respect to the conformal time 1], and primes mean 
derivatives with respect to r. 

As seen from the equations above, a Ist-order rotational perturbation cannot generate 
Ist-order perturbations of the mass-energy density and pressure. This is expected because 
p and P should be unchanged under the inversion of the rotation. For the same reason, 
iUi{r,ri) = lU^ir^rf) = 0. From Eq. (11), we see that if fi(r, r^) is independent of r, then 
iM^(r, 77) = fi(?7). That is, the fluid in the universe rotates with the metric at the same 



pace. Nevertheless, an r-dependent Q 



universe as that discussed in Ref. 



39 



r,p) allows us to discuss the centrifugal force for the 



401 ] for relativistic stars. Eq. (12) implies that Q{r,ri) 
must be in the form of a~^{r])A{r) + -8(77), where A{r) and -8(77) are arbitrary functions. 
Moreover, if the fluid is viscous, the R.H.S. of Eq. (12) will be equal to the Ist-order shear 
term of SnTro and this will free Q{r, rj) from the form above. 

Without loss of generality, we perform the following transformations: 

f{r,ri) = r'^a^{ri)Q'^{r,ri) — k{r,ri) — T{r,ri), (13) 
h{r,ri) = k{r,ri) — L{r,ri), (14) 

where T{r,r]) and L{r,ri) are arbitrary functions depending on f{r,r]) and h{r,r]). The 
first term of f{r,ri) comes from the transformation dO ^ dO — a{ri)Q{r,ri)dri. Using these 
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transformations to formulate the second-order EFEs, we find that 

2«2(r,?7) 



2Uo{r, Tj) 

2Ui{r, rj) 
2p{r, T}) 



2P{r, rj) 



2P{r,ri) 



2U3{r,r]) = 0, 
a{r]) 



I -2,U2n{r, T)) + ^i^^^ + k{r, rj) + T{r, rj) | , 



80^(77) - 40(77)0(77) 

+o(77)[-2^'(r,77)+L'(r,^)]}, 
1 



{-d(77)[2A;'(r,77) + 2T'(r,77)] 

4d(77)[3A;(r,77)-2L(r,77)] 



(15) 
(16) 

(17) 



{-4Aa^(77)[A;(r,r7) +T(r,77)] + 



327ra2(77) 
2[2A;'(r,77) -L'(r,77)] 



0(77) 



+ 12ra\'r])Q{r, 'r])Q'{r, 77) + r^a\'r])Q'\r, rj) - 4k"{r, rj) 



+2L"(r, 77) + Ar^a\ri)n(r, r])n"(r, 77) + 4 



Ao2(77) + 



o?{r]) 



-r'^a'^{r))Q,'^{r,r)) + 
1 



ra{r])fl{r, 77) — 



iU2{r,r]) 
rairj) 



}}, 



{k{T,r))+T{T,r)) 
(18) 



327rro4(77) ^^^^^'^^ [3a{ri)k{r, 77) - a{r])L{r, 77) + d{r))T{r, 77)] + 2o (77)[2rA;(r, 77) 
-rL{r, 77) + T'{r, 77)] + ro^(77)[-4AA;(r, 77) - 4AT(r, 77) + r^n'^{r, 77)]}, (19) 
^^^^^3^^^{-12d2(7/) [iU2{r, 1])^ + 4r2a^(77) [3a{r])k{r, 77) - a{r))L{r, 77) 
+a{rj)t{r, 77)] - r2o^(77)[4AA;(r, 77) + 4AT(r, 77) + 3rV{r, rj)] 
W(77){-4A[i«2(r, 77)]^ + 2r'[2k(r, 77) - L(r, 77) + T"{r, 77)]}}, (20) 
-{4ra(7/)d(7/)[3A;(r, 7;) — 2L{r,ri) + T(r, 77)] — ra^(7;)[4AA;(r, 7^) 



327rro4(77) 

+4Ar(r, 77) + r^n'^r, 77)] + 2a'^{ri)[2rk{r, rj) - 2rL{r, rj) + L'{r, rj) + T'{r, rj) 
+rL"{r, 77) + rT"(r, 77)]}. (21) 

Eqs. (19)-(21) are three different expressions for 2-P(?', ^) derived by the three 2nd- 
order spatial-spatial EFEs. Eq. (15) is expected by considering the symmetries of ^(odd)-, 
2;(even)-components of the four-velocity under the inversion of the rotation. 2U^{r,rj) — 
— 2^1 (r, 7/) 70^(77), which is non-zero in general and corresponds to the dynamical changes of 
2p(r, 77) and 2P{i^, v) r-dcpcndcnt rotational speed. In order to calculate these 2nd- 

order perturbations, we have to find the solutions of k{r,r]), L{r,ri) and T{r,rj). Since the 
pressure is the same along different directions at one point, Eqs. (19)-(21) are equivalent. 
Substracting Eqs. (19)-(21) from each other leads to two equations for solving -L(r, 77) and 
T(r, 77) when Q(r, 77) is specified while k{r, rj) is regarded as an arbitrary function indepen- 
dent of the rotation. The detailed derivations of these solutions are shown in APPENDIX 
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1. 



III. THE SACHS- WOLFE EFFECTS UP TO SECOND ORDER 

As the SW effect is invariant under the inversion of the rotation, the first non-zero SW 
effect due to rotation occurs in 2nd-order perturbations. The general formahsm of the 2nd- 



oder SW effect has 



the ideas in 



Q,[l5|, 



3een comprehensively discussed. In the following, we will make use of 



16 



171] and derive the 2nd-order SW effect of a rotating universe. 



The Cosmic Microwave Background (CMB) temperature observed at the origin towards 
a direction e can be written as 



Toie) 



-T,{p,d), 



(22) 



where u = —a~'^g^i,u^k^, the subscripts (O and e) denoting the origin and the last scattering 
hypersurface (LSH) respectively, = a^^{ri)dx^ / dr] is the four- velocity of the fluid in the 
comoving frame, k'^ = dx'^/dX is the wave vector of a light ray in the conformal metric with 
an affine parameter A, T^{p, d) is the temperature measured at the point p on the LSH, and 
d is the direction of the light (passing through the point p) observed at the origin. 

We show in Appendix II that the Ist-order SW effect due to rotation is zero and the 
2nd-order SW effect is 



Q\ sin 



2Ui{rx,'nx) . , , 2Uo{rx,'nx) 
- smd) + 



2A(i - 

1 r nla\r], 



a{ri\) a{r]x) 
[2d{ri,)T'{r„r],) - a{r],)L'{r„r]^)] 



k\X) 



\V0 

1% 



2 I 4A2(l-r]^)2 



T(-Asin0,A) , \ ■ a w ■ a ^(-Asin0,A) . 2 
h J (— Asm0, A)sm0 z sm 



d\ (23) 



2 V r, / r 2 

where = —77^ sin (j) and 7]^ denotes the conformal time of the last scattering. 

Eq. (23) determines the CMBA produced by the rotation of the universe once VL{r,7]) is 
specified. As an example, we consider the simplest case - stationary homogeneous rotation 
(i.e. f2(r, 77) = B{ri), B{ri) is an arbitrary function). Then, we have f{r,ri) = CQ^^^r"^ and 
h{r,ri) = where C is a constant. It is straight-forward to find that 

f = 0. (24) 



To explain this, we recall that i'U^(r, 77) = Q{r, 77) if Q{r, rf) is independent of r, which means 
that the fluid is rotating with the same phase as the metric. Therefore, the effect of the 
rotating metric cancels the relativistic Doppler effect caused by the sources rotating in a 
stationary metric. 

We make use of the previous example, i.e. f2(r, 77) = firotO^(^e)^^/(^e 0^(77)) with a = 
a?{ri^)/rl in Eqs. (A.10)-(A.18), to constrain the rotation of the universe. Using Eqs. (23), 
(A.IO), (A.16)-(A.18), we expand the CMBA as 
5T 

= 02 sin^ + 04 sin^ (p + Uq sin^ 

= AXi<l>, 0) + ^2K,°(0, 0) + A,Y,\(j), 6) + Ay6°(0, 9). (25) 

The values of A^s are listed in Table 1. 
Table 1 



-6.14188 
2 4.57532 
4 -1.67194 
6 0.25710 

We notice that the spherical harmonic expansion has non-zero coefficients only when 



771 = and even I for which Yi"^{(f),6) has cj. 
is unlikely to be related to the 'Axis of Evil' 



indrical and parity symmetries. The result 



y, y, 0, y, e 



lOj, a preferred direction of 



several low multipoles (especially quadrupole and octopole). However, in the general case, 
when we are located off the rotational axis, the cylindrical symmetry is broken and non-zero 
coefficients for other multipoles are allowed. For example, the CMBA on the two sides of 
the rotation axis will be affected differently by the rotation in general. Such an asymmetric 
effect enhances the dipole moment of the CMBA. Thus, its potential to explain the 'Axis of 
Evil' cannot be eliminated without further study. 

With = 0.742, Hq = 71.9km/s/Mpc and An ~ 10^^, we constrain Q^ot to be less than 
~ 6 X 10~^^ m in SI unit. That is, ^l{r^, rj^) is less than ~ 10~^ rad yr~^ in usual unit at the 
last scattering surface. Some CMBA maps generated with the rotation of the universe are 
shown in Fig. 1 as examples. Nevertheless, our result can be regarded as the first constraint 
of the rotation of a ACDM universe. 
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In Fig. 2, some normalized 2nd-order perturbed quantities along the light path of 
the last scattered photons are plotted as a function of ra{ri)Q{r, rj) with Q{r, rj) = 
Q^roto? {Tj^)r'^ / {rla? {rf)) and Vtrot ~ 6 x 10^^^ m. As expected, the perturbed quantities 
increase with the rotating speed. In Fig. 3, the angular velocity of matter iu'^{rri) and its 
difference from that of the metric along the light path of the last scattered photons are 
plotted against time. We can see that the angular velocity of matter can be negative while 
the rotation speed of the universe is always positive. Because of the r-dependence of fi(r, rj), 
the angular velocity of matter can be different from that of the metric in general as indicated 
in Eq. (11). These quantities are useful for studying the frame- dragging of the universe in 
the future. The distributions of the 2nd-order perturbed densities of matter are shown at 
two different times in Fig. 4. As shown in Fig. 2, 2W^(r, r/) is always positive, which means 
that matter is moving away from the rotating axis and hence the density is expected to be 
decreasing with time (shown in Fig. 4). 



IV. DISCUSSION AND CONCLUSIONS 



In this paper, we have developed a cosmological model that has a non-zero rotation but 
trivial shear in the Robertson- Walker framework with a cosmological constant. We have 
solved the EFE's up to 2nd-order perturbations of a flat ACDM universe with rotation as a 
Ist-order perturbation. We also set up the formulation for the 2nd-order SW effect due to 
the rotational perturbation and find that the effect only influences the spherical harmonics 
with even /'s. By making use of recent CMBA data, the angular speed of the rotation is 
constrained to be less than 10~^ rad yr~^ at the last scattering surface. 

The model of the universe here is different from the Bianchi models used in the literatures. 
First of all, our model is inhomogeneous with an axial symmetry in general while Bianchi 
models are homogeneous. Moreover, our model is shear-free and thus has the advantage 
that the SW effect and the constraint obtained are purely due to the global rotation. 

Conipared with previous works, the constraint is much weaker. For example. Barrow 
et al. 3| put a constraint of 1.5 x 10^^^ rad yr^^ on the rotaion of flat Bianchi models. 
This can be understood mainly because the effects of rotation in our model here show up 
as 2nd-order SW effects while in previous works they are Ist-order SW effects. For further 
study, we notice that our model here produces a 2nd-order outward radial velocity. It 
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may be used to explain parts of the accelerating expansion of the universe, and therefore, 
constraints on the global rotation can be obtained with Type la supernova data as proposed 
in 



34| . Although we only study a flat universe here, it is interesting to study the closed and 
open cases of our model in view of the significantly different constraints on closed and open 
Bianchi models j^, 27 1. 



The constraint obtained here, which is less than 10 ^ rad yr ^ at the last scattering 
surface, provides the first constraint of the shear-free rotation of a ACDM universe. 

APPENDIX I 

Since the pressure is the same along different directions at one point, Eqs. (19)- (21) are 
equivalent. By substracting Eqs. (19)-(21) from each other, we have 

T\r,i^)-rT"{r,r^) = 

+rKa^{ri)[?)Q! [r^r]) + ri7"(r, t])]^ + ra(?7)d^(r7){— 64r^a(?7)f2'^(r, 77) 
+3a(ry) [30'(r, r/) + r9."{r, r/)]^}}, (A.l) 
= 2rd^(?7)L(r, r]) + ra(r7)[— 4d(r7)L(r, 77) + 2d(?7)L(r, r])] 

+a^{r])[-2rM{r,r]) + r^fi'^(r, ?])] + 0? {r])[rL{r , 7]) - L'{r,r]) 
-rL"{r,r]) ~rT"{r,r])]. (A.2) 

These two equations are independent of k{r,ri). That is, k{r,ri) is only an arbitrary 
function unrelated to Q{r,ri) in general. Physically, k{r,r]) comes from the 2nd-order per- 
turbation of mass density (analogous to the diagonal perturbations of Schwarzschild metric). 
As we are interested in the effects of the rotation only, we set k{r,r]) = for simplicity. 

In this paper, we focus on a non- viscous fluid in the universe. Using the fact that 
op{ri) = A(l — Q/<^)/[S7TQAa^{ri)] for a flat ACDM universe and Eq. (8), we can simplify these 
two equations further 

T r, 7]) - rT r, r]) = — — — . 3. . — , A.3 

0^(7/) 2A{1 - n a) a^{r]) 

T d 

[a^(?7)L(r, r])] — a{ri)[L' [r^r]) +rL"(r, 77)] = ra{r])T"(r,r]) — {r]) A''^ (r) . (A. 4) 



a{'ri) drj 

Here, we remark that B{ri), which disappears from these two equations, is arbitrary and 
does not affect the Sachs- Wolfe effect. Such an arbitrariness is unrelated to the effects of 
the rotational universe on the CMBA. 
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Without loss of generality, we expand A{r) into Taylor series: A{r) — Qrot Yl'^=2 ^nf^ (we 
will explain below why it does not start from n — Qi). Substituting the series into Eq. (A. 3), 



T'(r,?7)-rT"(r,77) 

oo oo 

^ ^ /(n - / - VjCiCn-i-ir' 



2A(1 - n^)a^irj) 

where c„ = for n < 2. 

By separation of variables, we find that 



^-•j qo oo oo 

5^ + 2)(n - / + l)(n - / + 3)QC,_;+ir", (A.5) 



(A.6) 



where 



As Di{r]) and D2{i]) arc arbitrary functions of the homogeneous solutions for Eq. (A.5) and 
thus are independent of Q{r, rj), we are free to set them zero. The series of A[r) starts from 
n — 2 because the n — term can be absorbed into B{rj) while the n — 1 term, which 
produces ln(r) as the particular solution of T{r,r]), is rejected because of the singularity at 
r = 0. 

Similarly, we expand L{r,7]) = ^rotYl^=o ■^■n-{v)f^ substitute it into Eq. (A. 4): 

l{3n - l){n - I - 1) 



n=0 

oo oo 



{n-l)a%r)) ^'^""'"i^ 



n=5 1=2 

oo oo 



nljl + 2){n-l + l)in-l + 3)0a 
2A(l-QA)(n-l)a2(r?) 

where En{rj) = for n < 0. In general, by comparing the r" terms on both sides, we can 
obtain the recurrence relations for solving Enirf). As an example, we work out the simplest 
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case where il(r, 77) = ccQrot^^o ^{v) i^^ is a constant): 



= E2n+l{v), (A.ll) 

= ^^-^l^'ivWv)] - HV)E2{V): (A.12) 



ASriAa"^ 1 r 2 



a\7])E2{v)] - 160(77)^4(77), (A.13) 



A(l - l^A)a2(r7) a{7])dr]' 

[a'{ri)E,{ri)]-36a{ri)E,{ri), (A.14) 



14a2 1 d 



a^{r]) a{rj) dr] 

= ■^-^[a^(v)E2n-2m-'^n''a(v)E2n(v)- (A.15) 

We notice that there is a freedom to choose one of E2n{r]) arbitrarily, independent of 
Q(r, 77). To prevent infinite series , we set -£'6(77) = so that E2n{ri) = for ti > 3 and all 
non-zero E2n{'ny^ depend on (due to the rotation). We have 

MV) - -T7^dr,"dr,', (A.16) 



E2{V) 



° 1 



1 /-O 



- + 16a^(77")£;4(V') 



d77"d77', (A. 17) 



A(l - ^K)a{ri") 

Eo{v) = f^J\a\ri")E2{v")dri"dri', (A.18) 



which can be solved numerically. 



APPENDIX II 



Using Eq. (22) and expanding u^, k" up to 2nd-order in Vtrot as in Eqs. (4)-(7), we have the 
temperature anisotropy 

+[{ok^ lU^ + ik^ ou^) Uf + (oA:^ 2% + ik" + 2k'' ou^- (^-l) 

In order to calculate Eq. (B.l), we need to solve the geodesic equations for the light rays of 
the CMB, which are 

d ( dx''^ _ 1 dx'^dx'^ 
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As before, we expand the geodesic equations into different orders of Q(r, 77). To be consistent, 
gixv{x~') has to be expanded into 



For the zeroth-order, it is trivial that 

ok^ — (1, — sin 0,0, — COS0), 

where (j) is the zenith angle measured from the 2;-axis. 
For the first-order. 



(B.4) 



Ofl'22- 



dX 



+ 1920 



d_ 

dX 



dpX^' 

~d\ 

dix^ 
'~dX 



= 0, 
= 0, 



(B.5) 
(B.6) 



for i = 0,1,3. Therefore, ik^^{\) = {Aq, Ai,a{rix)^{rx,rix) — A2/r\,Az) where are 
constants to be determined. For simplicity, we assume that we are located on the rotating 
axis. Therefore, r]\ — X and r\ — —X sin (f) due to the cylindrical symmetry. Although the 
general case that we may be off the rotating axis is more realistic, the constraint here can be 
regarded as a good approximation provided that our distance to the rotating axis is small 
compared to that of the last scattering surface. The general case can be found by assigning 
a suitable dependence of 9 on rjx and fx- The term a(?7A)^(^A) ''7a) refers to the angular 
velocity of the comoving metric. 
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For the second-order, 

d ( dix^ 



— l5'20,0- 



d\X^ di^x^ 
dix^ d{^x^ ^ 2^700,0 ( dQX^\ ^ 2fi'ii,o / dt^x"^ 



d\ d\ 
d ( d2X^ 

dx V^''^ 



+ 



+ 2911 



2 V 

dox^ 



2^733,0 / doX' 



dX 



(B.7) 



dX 



09'22,l / dix 



2\ 2 



dX 



+ lfl'20,1" 



dix"^ d^x^ 2^7oo,l f dQX^\ 2fi'ii,i f doX 



dX dX 



dX 



1\ 2 



dX 



+ 



25'33,1 / 



3\ 2 



2 


\ 


d 




dX 


0922 


d 




dX 


09'33 



dX 



(B.8) 



-dx + ^^^^^ 

By setting k{r, rj) — and the requirement of a null geodesic {k^k^ = 0), we obtain 

^k^{X) = (2A;°(A),2A;^(A),0,0), 



where 



2k\X) 



-T'{rx',Vx')sin(f) + 



sin (f) 



2k\X) = 
To have a null geodesic 



T'{rx',Vx') . L'{ry,Vx') .2, f/ \ ■ , 
— ^-r H ^-T sm" (p - L[ry, r]y) sm 



dX', 
dX'. 



(B.IO) 

(B.ll) 
(B.12) 

(B.13) 
(B.14) 



= 22k\X) + T{rx,Vx) + '22k\X) sin (I) -L{rx, Tlx) siii" (I), (B.15) 

for A = 0. Prom the results of k^ and m^, we can easily verify the argument that the Ist-order 

perturbations of SW effect due to the rotation is zero and calculate the 2nd-order SW effect 
as 

ST, \ 

^ - • - sm0+ - - ■ - +2fc"(Aj + — 

Te I a{r]x) a{r]x) a{i]x) 

ft A sin ( 



2A(1 - Qa) 



aivx) a{r]x) 



0M^ , ^''-^'^^^>^) ,k\X) 



+ 



r(-Asin0,A) X • ^ u • ^ L(-Asin(/),A) . 2 
h i (— Asm0, A) sm0 sm ^ 



dA, (B.16) 
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where = 



—7]^ sin and denotes the occuring conformal time of the last scattering. 
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FIG. 1: The top map shows the effects of the rotation on bT^jT^ under the Mollweide projection 
with the z-direction pointing to (6, /) = (60°, 120°) with ^rot ~ 6 x 10^^^ m, which is the maximum 



ahowed by current CMB data. The middle map shows the original 5-year WMAP map 
the bottom map is a combined map of the two above. 



36( 1 ■ and 
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Log(i2(r,;7)//2(r,/7^)) 

FIG. 2: Some 2nd-order perturbed quantities are shown against ra{r))Cl{r,r)) along the hght paths 
of the last-scattered photons, where r\ = — Asin^ and r]\ = A. We set (f) = 7r/2. 
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FIG. 3: The velocity of matter iu'^{r\,'q\) and its difference from the rotation speed of the universe 
are plotted against time. Similar to Fig. 2, r\ = —A sine/), rjx = X and 4> = 7r/2. We note that 
the rotation of matter in the universe can be different and even opposite to the rotation of the 
universe. 



19 




1 1 1 1 1 1 1 1 1 
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FIG. 4: The distributions of the 2nd-order perturbed densities of matter are shown at two different 
time. The perturbed density is decreasing with r, which is consistent with the positive 2nd-order 
radial velocity {r, f]) ■ 
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